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We consider the dynamic critical behavior of the propagating mode for the order parameter 
fluctuation of the 0{N) Ginzburg-Landau theory, involving the canonical momentum as a degree 
of freedom. We reexamine the renormalization group analysis of the Langevin equation for the 
propagating mode. We Hnd the fixed point for the propagating mode as well as that for the diffusive 
one, the former of which is unstable to the latter. This indicates that the propagating mode becomes 
overdamped near the critical point. We thus can have a sufficient understanding of the phonon 
mode in the structural phase transition of solids. We also discuss the implication for the chiral 
phase transition. 
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The Quantum Chromodynamics (QCD) is believed to 
have a rich phase structure in the temperature (T) and 
baryon chemical potential (/tb) plane 0,13- The study 
of the dynamic (z.e., nonequilibrium) critical phenomena 
of the second order QCE)phase transitions including the 
chiral phase transition jj, II, IM la B B , the tricritical 
point, the critical end point iT^Tlini^ ll^ and the 
transitions associated with the color superconductivity 
[l^ is not only of theoretical interest but also of funda¬ 
mental importance for understanding of the relativistic 
heavy ion collision experiments and the early universe 
where the systems make time developments. 

One of the characteristic features of the dynamic criti¬ 
cal phenomena is the critical slowing down, which means 
that it takes very long time for the system near a critical 
point to relax into the equilibrium state. The long relax¬ 
ation time is attributed to the slow motion of the long- 
wavelength fluctuations of the so-called slow (or gross) 
variables. The slow variables are the fundamental de¬ 
grees of freedom for description of the dynamic proper¬ 
ties near critical points. Usually the slow variables are 
identihed with the collection of the order parameter and 
the conserved quantities of the system. 

The long-wavelength fluctuations of the slow variables 
make up what we call the slow modes. One should note 
that there are two kinds of slow modes; the propagating 
mode and the diffusive mode. The propagating mode 
involves oscillation with dissipation and corresponds to 
a pole with both the real and imaginary parts of the 
spectral function for the slow variables. The diffusive 
mode, on the other hand, is purely dissipative and may 
be called the relaxational mode, which corresponds to a 
pole with the imaginary part only. 

In general, it is possible to classify critical points into 
the universality classes. For the static (i.e., equilibrium) 
case, as is well known, the universality class can be de¬ 
termined solely by the symmetry and dimensions of the 
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system. For the dynamic case, a classification scheme was 
proposed by Hohenberg and Halperin M The dynamic 
universality classes are dependent on what kinds of slow 
variable (the order parameter and conserved quantities) 
are contained in the system as well as on the symmetry 
and dimensions. According to the scheme, Hohenberg 
and Halperin have classihed the whole critical points in 
the condensed matter physics in a lucid and systematic 
way [13. 

The dynamic universality class of the chiral phase tran¬ 
sition was first discussed in Ref. 0. Hohenberg and 
Halperin’s classification scheme tells us that the chiral 
phase transition belongs to the same dynamic universal¬ 
ity class as that of the antiferromagnet. In Refs. 0,0]: 
however, a crucial difference between the two systems has 
been pointed out. The difference is in the slow modes of 
the order parameter fluctuation. In the disordered phase, 
the order parameter fluctuation in the antiferromagnet is 
known to be a diffusive mode, while the order parameter 
fluctuation in the chiral phase transition gives the me¬ 
son or particle mode, which is a propagating mode. It 
has been argued (3 that in order to describe the meson 
mode appropriately, the canonical momentum conjugate 
to the order parameter in addition to the order param¬ 
eter itself and conserved quantities is needed as a slow 
variable. The necessity incorporating both the order pa¬ 
rameter and its canonical momentum may be accepted 
if we think of the Heisenberg equation for the Klein- 
Gordon field. We note that the canonical momentum 
is neither the order parameter nor a conserved quantity. 
This means that Hohenberg and Halperin’s classification 
scheme is not able to correctly classify the dynamic uni¬ 
versality class of the chiral phase transition and is not a 
complete one at its face value. 

In this Letter, we investigate the propagating mode for 
the order parameter fluctuation in the 0{N) Ginzburg- 
Landau theory that involves the canonical momentum as 
a degree of freedom. We analyze the Langevin equation 
for the propagating mode by using the renormalization 
group technique. For simplicity, we ignore other possible 
conserved quantities. The system we will consider is the 
simplest possible model. 
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In fact, we know a critical point in the condensed mat¬ 
ter physics, the soft mode for which is a propagating 
mode involving the canonical momentum. It is the struc¬ 
tural phase transition of solids M The soft mode is a 
phonon mode, which is a vibrational mode of the lattice. 
Experimentally, it is observed that the phonon mode be¬ 
comes overdamped, eventually turning into a diffusive 
mode near the critical point. This fact may suggest that 
the meson mode in the chiral transition also becomes 
overdamped and its dynamic universality class reduces 
to that of the antiferromagnet near the critical point as 
discussed in Ref. 0. This is indeed the case as we will 
see later. In a theoretical side, the phonon mode was 
examined in Ref. [H, in which essentially the Langevin 
equation for the propagating mode was analyzed. How¬ 
ever, the analysis will be found not to be adequate in 
the sense discussed below. Our investigation will eluci¬ 
date the profound background of the overdamping phe¬ 
nomenon. The purpose of this Letter is to give a reanal¬ 
ysis of the Langevin equation for the propagating mode 
of the 0{N) Ginzburg-Landau model and to examine the 
role of the canonical momentum in the critical dynamics. 
The model includes the structural phase transition as an 
N = 1 case and the chiral phase transition as = 4. 
Our analysis brings us new findings missed in the previ¬ 
ous work, which lead to a deeper insight into the struc¬ 
tural transition and the chiral transition as well as into 
the classification scheme of dynamic universality classes. 

The Langevin equation that describes the propagating 
mode for the order parameter fluctuation is given by a 
natural extension of that for the Brownian particle |l7| : 


d7ri(f, t) _ SF[(j)] _ d(j)^{x,t) 

dt 5(j)i{x,t) dt * ) ) 


{Ciix,t)Cjiy,t')) = kBT2TSij6{x - y)6(t - t'). 


( 1 ) 

( 2 ) 

( 3 ) 


The (l)i{x, t) is the A^-component (i = I, 2, • • • ,N) order 
parameter and {x, t) is its canonical momentum. The 
F[(j)] is the usual 0{N) symmetric Ginzburg-Landau free 
energy given by 



y(j)i{x,t) ■ V(j)i{x,t) -l-ro(/>^ -I- 


( 4 ) 

where d is the spatial dimension and rg and u are the 
usual static parameters. The dynamic parameters A and 
r represent the square of the propagating velocity and 
the damping constant, respectively. The C,i{x,t) is the 
white noise which satisfies Eq. ®. 


In the Fourier space, the Langevin equation reads 
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where is the system volume. The wavenumber of the 
fluctuations is cutoff at A. 

Our Langevin equation for the propagating mode 
should be in some connection with the Langevin equa¬ 
tion for the diffusive mode. In general, when the fric¬ 
tion or damping constant is very large, the oscillation 
becomes overdamped and the oscillatory nature is lost. 
The Langevin equation for the propagating mode should 
be reduced to that for the overdamped or diffusive mode 
for the large damping constant. Actually, this reduction 
can be proven explicitly if the nonlinear interaction cou¬ 
pling u is absent ^3. Consider the Langevin equation 
(IHTtITIi with u = 0. If we take the overdamped limit by 
imposing the condition AT^ ^ rg + k^, then the canonical 
momentum turns out to be the faster degree of freedom 
and we can integrate it out explicitly to find the Langevin 
equation for the slower degree of freedom 

^ k<A 

(^')> = kBT2jS,,6j: _j^,S{t - t'), (10) 

where is the renormalized noise term and 7 = l/T 

1^. This equation is for the diffusive mode and gives 
nothing but the model A with u = 0 in Ref. M 

Now we apply the renormalization group to the 
Langevin equation for the propagating mode, not neces¬ 
sarily assuming the overdamping condition. The renor¬ 
malization group analysis leading to the recursion re¬ 
lation 11511811 has alreadv been oerformed in Ref. M 
Although the following calculation is not new, we will 
present the calculation below to make the discussion self- 
contained. 

The renormalization group program consists of the two 
procedures mill; (i) integrate out the short-wavelength 
fluctuation with A/b < k < A, and (ii) make the scale 
transformation of After the pro¬ 

cedures, we have the recursion relation for the parameters 
(A, T, rg,u). We employ the e-expansion assuming that u 
is of order e, where e = 4 — d. In this Letter, we consider 
only the lowest order of the expansion. The dynamic 
response function for the order parameter fluctuation is 
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given in the form of 


Gik^uj) 


1 

— + ro + — iwr + S(fc, w) 


( 11 ) 


where Il{k,uj) is the IPI self-energy. 

The diagrammatic rules for the perturbation are al¬ 
most the same as those for the diffusive model (model 
A in Ref. ^3)- Only Ih® difference is the form of the 
response function. For the propagating mode, it includes 
the term as seen from Eq. CD- For the overdamped 
limit where AF ^ w (or -s/AF ^ A), the term disap¬ 
pears and the response function reduces to that of the 
diffusive mode We note that for the particle mode, 
which is a propagating mode and described by the system 
with A = 1, the pole position lies in the time-like region 
in the u-k plane. This is contrary to the diffusive mode 
for which the pole is always in the space-like region. 

After the integration of the fluctuation, the new pa¬ 
rameters are defined by 


Tq = lim G ^{k,uj), 

k ^ u )—»-0 


F' = 
1 

A “ 


d 


lim . 

k,u/^0 0( — lUJ) 
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G-^{k,uj), 
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( 12 ) 

(13) 
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where the last equation is absent in the diffusive model 
and new in the present model. Up to the lowest order 
in the e-expansion, the self-energy is given by the tad¬ 
pole diagram, which has no frequency and wavenumber 
dependence. Thus F' and receive no corrections from 
the fluctuation integration. After the scale transforma¬ 
tion, we obtain the recursion relation; 

V = ‘'"“v 

r' = (16) 

r'o = 6^(ro -k Ar), (17) 

u' = b^-^-^(u + Au). (18) 

We note that Eq. CD is given in Eq. (4.36) in Ref. 0 

The first two equations are for the dynamic parameters 
while the last two for the static ones. We note that 
the recursions for the dynamic and static parameters are 
decoupled to this order. The recursions for the static 
parameters are the usual ones for the static Ginzburg- 

Landau theory with Ar = u{N/2 -k 1) / ( 27 rU 
and Au = —u^(A-k8)/2 f jj^(ro + and give rise 

to the Gaussian and Wilson-Fisher (WF) fixed points as 
usual. Although there are two dynamic parameters, one 
of them just fixes the time scale and only some ratio of 
the two is meaningful. It is useful to define the dynamic 
parameter as p = -s/AF/A. From the recursion relation, 
we find that there arise two fixed points associated with 



FIG. 1: The two fixed points and the renormalization group 
flow. 


the WF fixed point in the parameter space (o. rn. 

See Fig. Q] 

The two fixed points are located at p = 0 and p = oo, 
which will be denoted as P and D, respectively. The fixed 
point P has the dynamic critical exponent z = 1 and 
represents the purely propagating mode, while the fixed 
point D has the dynamic exponent z = 2 pn| and cor¬ 
responds to the diffusive mode. The dynamic exponents 
are obtained in a usual manner |l8ll2(l| . The arrows indi¬ 
cate the flow of the renormalization group. We see that 
the parameter p is relevant with respect to the fixed point 
P and the line extended from the Gaussian fixed point 
to the WF fixed point gives the critical ‘surface.’ Thus 
the fixed point P is unstable to the fixed point D. The 
system represented by the point near the critical surface, 
such as the cross in Fig. ^ is firstly taken to the vicin¬ 
ity of P but eventually driven away to the fixed point D; 
namely, there occurs a crossover between the two fixed 
points. 

An illuminative example of a crossover phenomenon 
between two fixed points is provided us with the static 
critical phenomena in the magnetic systems 0- Gon- 
sider the Heisenberg ferromagnet which has the rota¬ 
tional 0(3) symmetry. Its critical property would be 
controlled by the associated fixed point, i.e., the Heisen¬ 
berg fixed point. The full rotational symmetry may be 
broken into the uniaxial one by, say, distortion of the lat¬ 
tice, the effect of which is described by the appropriate 
anisotropic Hamiltonian with the coupling constant g. It 
is known that the Heisenberg fixed point is unstable 
with respect to the parameter g. The stable fixed points 
are supplied by the Ising or XY fixed points, depending 
on the sign of g. A point with small g in the parameter 
space is driven by the renormalization group to the Ising 
or XY fixed point via the vicinity of the Heisenberg fixed 
point. This means that the system firstly exhibits the 
critical behavior of the Heisenberg universality class, but 
the Ising- or XY-like behavior eventually shows up as the 
critical point is approached. 

Our analysis indicates that the propagating mode 
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crosses over into the diffusive mode near the critical 
point: At hrst, the propagating mode is softened with 
the exponent z = 1 under the control of the propagat¬ 
ing hxed point, and then it becomes overdamped to be 
governed by the diffusive fixed point with z = 2. 

Now we have obtained a lucid explanation for the 
phonon mode in the structural phase transition : The 
phonon mode changes its behavior from the propagating 
into the diffusive as a consequence of the crossover be¬ 
tween the two fixed points. As we mentioned, the re¬ 
cursion relation 11511811 was alreadv derived in Ref. 0 
However it was utilized only to discuss the stability of 
the diffusive hxed point with respect to the small per¬ 
turbation in the p direction. Namely only the vicinity of 
the diffusive hxed point was investigated. In the present 
analysis, we have examined the whole region in the pa¬ 
rameter space within the validity of the e-expansion, and 
found out the unstable propagating hxed point. 

The hnding of the unstable propagating hxed point 
brings us new insights into the structural phase transi¬ 
tion: We can have a simple understanding of the over¬ 
damping phenomenon as a crossover between the two 
hxed points. The mechanism of how the propagating 
behavior changes into the dihusive one has become more 
concrete. It is clear that because of the instability of the 
hxed point, the propagating mode inevitably becomes 
overdamped even if we start from anywhere in the pa¬ 
rameter space except in the p = 0 plane. Moreover the 
universality of the propagating behavior of the phonon 
mode are conhrmed by the hxed point. These ingredients 
give a full explanation in terms of the renormalization 
group language, which is sufficient for understanding of 
the phonon mode in the structural phase transition |22l| . 

We note that it is not until the propagating hxed point 
is given that its universality is guaranteed. In fact, the 
renormalization group analysis of the phonon mode has 
been restricted only to the diffusive behavior near the 
critical point up to the present. The propagating hxed 
point we have found gives us a hrm basis to discuss 
the universal nature of the propagating behavior of the 
phonon mode. 

It should be noted that though being unstable, the 
propagating hxed point gives one dynamic universality 
class in the sense that hxed points and universality classes 
make a one-to-one correspondence. The Heisenberg hxed 
point in the magnetic system actually gives one univer¬ 
sality class although it is unstable against the anisotropic 
interaction. As can be seen from the recursion relation 
JTM . within the p = 0 plane including the propagating 
hxed point, the renormalization group how closes itself. 
If the system exists which has the parameters just in the 
plane, the propagating hxed point is no longer unstable 
. If the system takes the parameters very close to the 
plane initially, the overdamped region near the critical 
point is correspondingly small and the critical behavior 
is almost governed by the propagating hxed point. Thus 
the propagating behavior constitutes one dynamic uni¬ 
versality class. 


Since the propagating mode eventually reduces to the 
diffusive mode near the critical point, the canonical mo¬ 
mentum becomes a rapid degree of freedom and fades 
out of the member of the slow variables at that point. 
Thus if one restricts one’s interest to the very vicinity 
of the critical point, Hohenberg and Halperin’s classih- 
cation scheme of dynamic universality classes, which is 
based on the order parameter and conserved quantities, 
practically works. However, the canonical momentum is 
still an important degree of freedom to hnd the dynamic 
universality class associated with the propagating hxed 
point. In this respect, Hohenberg and Halperin’s classi- 
hcation scheme may be regarded as incomplete. 


We now turn to the chiral phase transition. Our anal¬ 
ysis predicts that as the system approaches the critical 
point, the meson mode turns into the diffusive mode after 
the softening [^. Thus the dynamic universality class 
of the chiral transition certainly reduces to that of the 
antiferromagnet as argued in Ref. Q . However the prop¬ 
agating behavior of the meson mode is still important to 
analyze since it shows universal properties belonging to a 
hxed point. There have been various approaches, such as 
the mode coupling theory 0,0, the Nambu-Jona-Lasinio 
model 0 and the microscopic approach in Ref. 0 . Un¬ 
fortunately, the overdamping phenomenon was not no¬ 
ticed in these analyses within the adopted approxima¬ 
tions. Actually, up to now there have been two kinds of 
works based on different standpoints: Some works treat 
the meson mode as a diffusive mode while others assume 
the mode to be propagating. These two standpoints are 
reconciled now that the crossover between the two kinds 
of modes has been realized. We should note that the 
overdamping of the meson mode itself has a significant 
physical meaning that the sigma mesons and pions are 
not able to propagate and lose a particle nature near the 
chiral phase transition. 


To summarize, we have found the fixed point for the 
propagating mode, which is unstable to the fixed point 
for the diffusive mode. This means that the propagating 
mode for the 0{N) Ginzburg-Landau theory becomes 
overdamped near the critical point. The analysis gives 
a satisfactory account of the character change of the 
phonon mode in the structural phase transition and also 
predicts the fate of the meson mode near the chiral phase 
transition. In the future work, we will investigate the 
higher order calculation in the e-expansion, which will 
clarify the deviations of the dynamic critical exponents 
from the mean held values obtained in this work. The 
1/A^ expansion is also interesting. Moreover it is neces¬ 
sary to take account of the conserved quantities for the 
proper description of the structural phase transition, the 
chiral phase transition and so on. We hope to report the 
progresses in these directions in the future. 
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